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ABSTRACT 

The Cygnus region is a very bright and complex portion of the TeV sky, host to unidentified 
sources and a diffuse excess with respect to conventional cosmic-ray propagation models. Two of 
the brightest TeV sources, MGRO J2019+37 and MGRO J2031+41, are analyzed using Milagro 
data with a new technique, and their emission is tested under two different spectral assumptions: 
a power law and a power law with an exponential cutoff. The new analysis technique is based 
on an energy estimator that uses the fraction of photomultiplier tubes in the observatory that 
detect the extensive air shower. The photon spectrum is measured in the range 1 to 200 TeV 
using the last 3 years of Milagro data (2005-2008), with the detector in its final configuration. 
MGRO J2019-t-37 is detected with a significance of 12.3 standard deviations (a), and is better 
fit by a power law with an exponential cutoff than by a simple power law, with a probability 
> 98% (F-test). The best-fitting parameters for the power law with exponential cutoff model 
are a normalization at 10 TeV of 7^1 x 10"^" s"^ m"^ TeV"\ a spectral index of 2.0ti'^ and 
a cutoff energy of 29^^^ TeV. MGRO J2031-h41 is detected with a significance of 7.3a, with 
no evidence of a cutoff. The best-fitting parameters for a power law are a normalization of 
2.4^0,1 X 10-^° s^i m-2 TeV^i and a spectral index of 3.08j:[J The overall flux is subject 
to an ~30% systematic uncertainty. The systematic uncertainty on the power law indices is 
~0.1. A comparison with previous results from TeV J2032-f4130, MGRO J2031-h41 and MGRO 
J2019-I-37 is also presented. 

Subject headings: Gamma rays: observations 
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1. Introduction 

The Cygnus region is a part of our Galaxy 
of active massive star formation and destruction, 
and has been studied over a broad range of wave- 
lengths, including radio, X-ray, GeV and TeV 
gamma-ray, as well as in cosmic rays. From GeV 
up to multi-TeV energies the Cygnus region is the 
brightest diffuse gamma-ray source in the north- 
ern hemisphere ( Hunter et al.lll997h . 

One of the challenges in analyzing the Cygnus 
region at TeV energies, is the proper separation 
of the gamma-ray flux that is attributed to the 
point or extended sources in the region or to the 
diffuse emission. Previous Milagro analyses com- 
puted the diffuse emission from the region using 
two alternative methods to isolate the contribu- 



Cygnus region (I € [65°, 85°] and b G [-2°, +2°]) 



(lAbdo et al.ll2nn9h 



tion from the resolved sources (jAbdo et al.l 12007 . 
2008I ). and found that at TeV energies the flux is 
still in excess with respect t o the predict ed flux 
from the GALPROP model (IGALPROPI) . Mila- 
gro also published the discovery of two unidenti- 
fied TeV sources in the region, MGRO J2 19-K3 7 
and MGRO J2031+41 (|Abdo et al.ll2007al l2009^ . 
The location of MGRO J2019-h37 was found to 
be consistent with two EGRET sources, while the 
best fit position for MGRO J2031+41 was near 
two EGRET sources and the HEGRA unidentified 
source TeV J2032-I-4130. In a correlation study 
connecting the TeV sources discovered by Mila- 
gro with sources detected above lOcr (the so-called 
Bright Source List, BSL) by the Fermi Large Area 
Telescope (LAT), the two aforementioned bright- 
est Milagro sources in the region were associated 
with GeV pulsars (|Abdo et al.ll2009cl l?h. MGRO 
J2019-h37 and MGRO J2031-h41 are currently as- 
sociated with two and one pulsars identified b y the 
Fermi LAT, respectively (|Abdo et al.ll2011ct ). Mi- 
lagro also detected two candidate sources in the 
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Recently, several collaborations have presented 
new surveys of the Cygnus region, confirming 
the complexity of the region and showing the 
highly structured diffuse emission. The TeV 
emission from the position of MGRO J2019+37 
has b een confirmed b y the VERITAS experi- 
ment llAliu et al-lboill ). The VERITAS spatial 



counts map shows a clear structure associated 
with MGRO J2019-H37, a region of extended 
emission which seems to be produced by previ- 
ously unresolved sources. At lower energies, the 
Fermi LAT collaboration has also recently pub- 
lished new results on the diffuse emission from the 
Cygnus region. An extended source, the so-called 
Cocoon, overlapping with MGRO J2031-I-41, has 
been found and its emission has been explained by 
freshly accelerated cosmic rays, trapped in a shell 
of photon-dominated emission formed by stellar 
winds and su pernovae, as shown by mid- infrared 
observations ( Abdo et al. 2011bl ). The spectrum 
from this region is hard, extending up to 100 GeV 
with no evidence of softening, and could be ex- 
plained as the result of pulsar-accelerated particles 
or as an active super-bubble. The average diffuse 
emission from the region has also been ana lyzed 
at MeV to GeV energies (|Abdo et al.ll2011al ). De- 
spite the very rich source population, this emission 
is similar to that of the local interstellar space, and 
there is no evidence that it is necessary to include 
an extra contribution in the model, resembling 
the diffuse excess previously measured by Mila- 
gro (|Abdo et al.ll2007l \200& }. Most recently, the 
ARGO-YBJ collaboration presented the results 
of a survey of the Cygnus region in the energy 
range of 600 GeV to 7 TeV. MGRO J2031+41 
is observed with a significance of 6Aa and a flux 
consistent with previous Milagro results, but there 
i s no significant de tection of MGRO J2019+37 
(|Bartoli et al.ll2012[ ). 

Here, we present a new analysis of the last 
three years of data collected with the Milagro 
experiment (2005-2008). An improved gamma- 
hadron separation and a newly-developed tech- 
nique are exploited to reconstruct the energy 
spectra of gamma rays from the sources in the 
Cygnus region. The best fits for the spectra of 
the two brightest Milagro sources are presented: 
MGRO J20194-37, a source with a post-trials sig- 
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nificance in excess of 12cr between 1 and 100 TeV, 
and MGRO J2031+41 a sour ce with a p ost-trials 
significance in excess of la (jAbdo et a l. 2007a). 
We compare our spectra with results from the 
HEGRA, M AGIC, Whippl e and AR GO-YBJ ex 



periments dAharonian etTal . 200 5^: jAlbert et al 
20081 : iLang et al.ll2004HBartoli et allf2012[ ). 



2. Analysis Technique 

The Milagro detector was located in the Jemez 
Mountains in New Mexico, at an altitude of 2,630 
m a.l.s. It was operated from 2001 to 2008, and, 
in its final configuration, consisted of two com- 
ponents: (1) a central pond (60 x 80 m^, 8 m 
deep) with two photomultiplier tube (PMT) lay- 
ers, a shallow air shower layer consisting of 450 
PMTs, and a deep muon layer consisting of 273 
PMTs and (2) an array of 175 single-PMT tank 
detectors surr ounding the poud cover ing an area 
of 40,000 m2 (jAtkins et al.ll200ll |2003| ). 

A detailed description of the analysis method 
and the parameters used here can be found in 
the Milagro paper on the spectral m easurement of 
the Crab Nebula (jAbdo et al.ll201lh . Systematic 
uncertainties of this method are estimated to be 
~30% on the overall flux and ~0.1 on the spectral 
index. 

The new element in this an alysis, with re- 
spect to the previous approach (.Abdo et al.ll2007al 
20091) . is the introduction of an estimator, for 
the energy of the primary particle initiating the ex- 
tensive air shower, based on the number of PMTs 
hit for each event. The J- parameter is defined as 
the sum of two fractions: 
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where Na.s. and Nt.a. are the number of PMTs 
detecting the event, while Nf;i and NJf^X. 

are the 

number of functional PMTs in the air shower layer 
and in the tank array, respectively. Because the 
typical energy resolution of Milagro is quite broad 
(50-100%), and the median energy associated with 
a given J- value is dependent on the spectral as- 
sumption for the observed source, the fit is per- 
formed in the measured !F space rather than in 
the energy space. 

The reconstructed Milagro data contain infor- 
mation about the direction and J- value of air 



shower events. According to their F value, events 
are sorted into 9 bins (0.2 < J" < 2), result- 
ing in a total of 9 signal skymaps, binned in 
0.1° X 0.1° pixels. For each F bin, background 
maps are calculated using the Direct Integration 
method (jAbdo et all l201l[) . The two brightest 
regions of interest in the sky, a 2° x 2° area 
around the position of the Crab Nebula and a 
band around the Galactic plane (—2.5° < b < 
2.5°) ar e excluded when calculating the back- 
ground (jAbdo et al.ll201l[ ). Rather than discrim- 
inating between gamma-ray and cosmic-ray initi- 
ated air showers, a weight is applied to all mea- 
sured events, in both signal and background maps, 
where gamma-ray like events are assig ned a higher 



weigh t than cosmic-ray like events (jAbdo et al 
20111 ). The gamma-ray excess with respect to 



the estimated background is calculated as the dif- 
ference between signal weights and background 
weights, and is filled into the so-called excess map. 
Based on these excess maps we compute the F dis- 
tributions used in the energy fit. 

Excess maps are smoothed according to the de- 
tector angular resolution (or point spread func- 
tion, PSF), which is a function of F and varies 
between 0.3° and 0.7°. The measured F distribu- 
tion for the target source in a 0.1° x 0.1° bin, as a 
result, is the average excess coming from a PSF- 
wide region around the nominal source position. 
Since the background is measured with the direct 
integration method removing the Galactic plane 
{b > 2.5° and b < —2.5°), the measured excess 
from sources in the Galactic plane includes any 
other diffuse or extended emission possibly present 
in the vicinity of the source. The Galactic diffuse 
background was estimated to contribute up to 15% 
of the fl ux at 35 TeV for t he weakest Galactic BSL 
source (jAbdo et al.ll2009l) . 

We test the emission from the two sources for 
two different spectral hypotheses: a power law, 
defined as 
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where iVioTeV is the normalization scale set at 10 
TeV and a the spectral index, and a power law 
with an exponential cutoff, defined as 



dN _ E 
d^-^^"^'=^llOTW 



exp(-^), (3) 
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where Ec is the cutoff energy. 

The fit to the data is performed comparing the 
measured excess to simulations. First, a set of sim- 
ulated data is generated varying sensitive parame- 
ters (i.e. A^ioTcVj and Ec). The best spectral pa- 
rameters and the corresponding fit probability are 
then found using a minimization, comparing 
the measured and the simulated J- distributions. 
Uncertainties on fit parameters are computed us- 
ing 1-CT contours of histograms (as discussed in 
sections l3.ll and 13. 2p . The uncertainty is defined 
as the distance between the best fit value and the 
lower and upper edges of the 1-a contour. 

3. Results and Discussion 

The results presented here are obtained from 
the analysis of the Milagro data taken with the 
detector in its final configuration. The start date 
is 2005-10-22 20:39:16 GMT and the stop date is 
2008-04-15 00:02:53 GMT, corresponding to a to- 
tal of 906 days of observation, resulting in 832 in- 
tegrated days after data quality cuts. 

The positions of MGRO J2019-I-37 and MGRO 
J2031-f41 are obtained fitting Milagro excess 
maps with a 2-dimensional Gaussian function with 
equal widths in right ascension and declination 
i'^RA — <^Dec), plus a Constant base to take into 
account contribution from diffuse emission. The 
fit is performed in a square region of ±3° around 
the excess. Table [1] shows the positions used in 
this paper and a comparison to other surveys. 

3.1. MGRO J2019+37 spectrum 

In Fig. [1] the distribution of data is shown 
with the two simulated J- distributions for the best 
fitting power law and power law with exponen- 
tial cutoff models. These distributions show the 
weighted excess coming from the source position 
after background subtraction. The excess from 
MGRO J2019-I-37 is compared to the excess mea- 
sured from the entire Cygnus region (65° < / < 
85° and -2° < 6 < 2°). We find that the emission 
from MGRO J2019-H37 accounts for -12% of the 
excess from the entire region. 

Figs. [2]and[3]show the 1 and 2-a contours from 
the fit for the two tested hypotheses. 1-a con- 
tours are used to compute uncertainty on the pa- 
rameters as described in section [2] 



The best fit spectra for the two tested hypothe- 
ses (Eqs. [2]and|3l) are shown in Fig. SI The best-fit 
spectral parameters are summarized in Table [2] 

With the power law hypothesis, the best-fit pa- 
rameters for the normalization and the spectral in- 
dex are iVioTcV = (4.1±0.6)-10-i° s"! m'^ TeY~^ 
and a = 2.78 ± 0.10 respectively. The is 16.12 
for 7 degrees of freedom (d.o.f.), which gives a x^ 
fit probability of 0.024. 

In the case of the power law with cutoff hy- 
pothesis, the best fit parameters are iVioxcV — 
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10 



-10 



29l^^ TeV, 



s-i m-2 TeV-i, a = 2.0;?;^ and E^ = 
The x^ is 1-924 (6 d.o.f.), which gives 



a X fit probability of 0.93. 

The x^ fit probabilities suggest that the power 
law with a cutoff model (x^ probability — 93%) 
fits the data better than a simple power law 
(x^ probability = 2.4%) W e use the F-test 
( Bevington and Robinsonll2003[ ) to compute if the 
improvement in the fit is significant. For MGRO 
J2019-I-37, the simple power law is disfavoured at 
a C.L. > 98%. 

Previous Milagro analyses quoted the fiux of 
MGRO J2019-H37 at 20 and 35 TeV respectively 
(|Abdo et al.ll2007al lioogh . Those values, using a 
different analysis technique, are in agreement with 
the r esults presen ted here. An independent anal- 
ysis ( AllenI 20071) . which used a different energy 
estimator for the particle initiating the air shower 
and a different parameter to distinguish between 
gamma and cosmic rays, _ also is consistent with 



our results. ARGO-YBJ ( Bartoh et al.ll2012D has 



not observed any significant emission from MGRO 
J2019+37 in the range 600 GeV to 7 TeV. The ex- 
planations put forward in (jBartoli et al.ll2012l ) are 
an insufficient exposure above 5 TeV or a possible 
time variability of MGRO J2019-I-37. The 90% 
C.L. upper limit from ARGO-YBJ is consistent 
with our best fitting model, a power law with an 
exponential cutoff (see Fig. |4]). On the other hand, 
the simple power law model (a = 2.78 ± 0.10), al- 
ready disfavored by the F-test, does not agree with 
the ARGO-YBJ resuffs. 

3.2. MGRO J2031+41 spectrum 

The J- distribution from data and those sim- 
ulated for the two best fit models are shown in 
Fig. m MGRO J2031-I-41 accounts for -6% of 
the excess from the entire Cygnus region. 
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Fig. 3.— MGRO J2019+37: 1 and 2-a C.L. con- 
tours for the power law with cutoff fit. left 'panel 
- NioTeV vs a projection, right panel - Ec vs a 
projection. The points indicate the best fit result 
(minimum x^)- The increments for the 1 and 
2-cr contours are 3.53 and 8.02 respectively. 



Fig. 1.— MGRO J2019+37: distributions of the 
parameter used to estimate the photon energy, T. 
The spectrum fit is performed in the T space. 
Black squares represent the data. Red triangles 
and blue circles represent the simulated distribu- 
tions for the best fit assuming a power law with a 
cutoff and a simple power law, respectively. The 
unit on the y axis is the weighted excess per day, 
the unit on the x axis is the value. 
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Fig. 2.— MGRO J2019-t-37: 1 and 2-cr C.L. con- 
tours for the power law fit. The point indicates 
the best fit result (minimum x^). The incre- 
ments for the 1 and 2-cr contours are 2.30 and 6.18 
respectively. 
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Fig. 4.— MGRO J2019-1-37: energy spectra. The 
best fit is obtained for a power law with cutoff 
model (in red). The power law model is also shown 
(in blue). The shadowed area represents the f-cr 
band, obtained by varying the parameters within 
the 1-CT contour. ARGO-YBJ 90% C.L. upper 
limits for MGRO J20f9+37 are shown in black 



(jBartoh et al.ll2012l) 
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Fig. |6] and [7] show the 1 and 2-(t contours from 
the fit for the two tested hypotheses. 

Best-fit spectral parameters are summarized in 
Tabled 

For the power law hypothesis, the best fit pa- 
rameters are TViotcV = 2.41° j-10-i° s'^ m'^ TcV" 
and a = S.OStg j?. The is 6.98 (7 d.o.f.), with 
a fit probability of 0.43. 

For the power law with cutoff hypothesis, the 



+90.0 
-1.9 



10 



-10 



44t:^ TeV. 



best fit parameters are A^iotcV = 3.7 
s~i TeV"\ a 
In this case we are not able to constrain the upper 
limit of the 1-cr uncertainty for the cutoff energy 
(see right panel of Fig. [7]). The maximum Ec 
simulated value is 1000 TeV, almost an order of 
magnitude above the highest detected energy. The 
X^ is 5.80 (6 d.o.f.), with a x^ fit probability of 
0.45. The best fit spectra for the two models are 
shown in Fig. [8] 

The improvement in the fit obtained with the 
more complex model is not significant, according 
to the F-test, and the simple power law is to be 
preferred over the power law with the exponential 
cutoff. 

The flux of MGRO J203H-41 at 20 and 35 
TeV, measured by pre vious Milagro analyses 
(|Abdo et al.l2007all2009l ) , is in agreement with the 
results presented here. Results from ARGO-YBJ 
MAGIC, HEGRA and Whipple ( | Bartoli et al 



2012t lAlbert et al.l l2008t lAharonian et al 



2005t 



Lang et al.l |2004[ ) are shown in Fig. El The 
HEGRA and MAGIC models are mutually con- 
sistent, but they disagree with our best fit in 
terms of integral flux in the overlapping energy 
range. Between 1 and 10 TeV the flux measured 
by MAGIC accounts for only '^3% of the fiux 
measured by Milagro. The spectrum as measured 
by the air Cherenkov telescopes (ACTs) is much 
harder (a ~2) than the Milagro power law best 
fit (a '^3). This discrepancy can be explained by 
the following two facts. First, the angular reso- 
lution of HEGRA and MAGIC (< 0.1°) is much 
better than that of Milagro (0.3° to 0.7°). The 
Whipple flux was measured at 0.6 TeV, with a 
PSF of 0.21°, and it agrees better with the ex- 
trapolation of the Milagro result at lower ener- 
gies. The same is true for ARGO-YBJ, which, 
with an angular resolution between 0.47° and 
2.8°, measured an emission consistent with the 



results presented here. The measurement pre- 
sented in this paper accounts for photons coming 
from a larger region around the nominal source 
position compared to ACTs. Second, the way the 
background is computed is also different. Whip- 
ple used the ON/OFF method and observed the 
source with a significance of 3.3cr, while MAGIC 
primarily was operated in wobble mode (5.6cr). 
A possible explanation of the presented results is 
that TeV J2032-I-4130, whose extension is slightly 
larger than the HEGRA and MAGIC angular res- 
olution, is surrounded by an extended emission, 
^4° X 3° according to the Milagro map (see Fig. |9] 
discussed in Sec. 13. 3|) . Therefore, Milagro, Whip- 
ple and ARGO-YBJ are not able to disentangle 
the extended emission from the central source and 
observe a higher flux. MAGIC and HEGRA, on 
the other hand, take into account the extended 
emission as a background. As a result, the ACTs 
measure a fainter emission and a harder spectrum 
for TeV J2032-I-4130. 

3.3. Morphology 

The significance map of the Cygnus region for 
the last three years of Mflagro data (2005-2008) 
is shown in Fig. [9l An energy-dependent PSF 
smoothing is applied. The top panel shows the en- 
tire region in Galactic coordinates. Milagro posi- 
tions for MGRO J2019+37 and MGRO J2031+41 
are shown, as well as the high-significance Fermi 
LAT sources (> lOcr) and the extended Fermi Co- 
coon. 

The Fermi counts map for photon energies 
above 10 GeV overlaid with the Milagro signif- 
icance contours in the Cygnus region is shown in 
Fig. [TOl Milagro detects an excess consistent with 
the position of the two pulsars PSR J2021-I-3651 
and PSR J2032-h4127, associated with MGRO 
J2019-h37 and J2031-H41 respectively, but no de- 
tection of significant emission coincident with PSR 
J202H-4026. 

MGRO J20 19-1-37 is observed with a signifi- 
cance of 12.3(7. Its extension is obtained with the 
2-dimensional Gaussian fit discussed in Sec. [3l 
The fit results in cr = 0.75°, slightly larger then the 
angular resolution of the detector in its final con- 
figuration. There is evidence from VERITAS that 
the extended source Milagro observes is the result 
of the superimpo sition of close point-like sources 
(|Aliu et al.ll20li . 
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Fig. 5.— MGRO J2031+41: distributions of the 
parameter used to estimate the photon energy, T. 
The spectrum fit is performed in the T space. 
Black squares represent the data. Red triangles 
and blue circles represent the simulated distribu- 
tions for the best fit assuming a power law with a 
cutoff and a simple power law, respectively. The 
unit on the y axis is the weighted excess per day, 
the unit on the x axis is the T value. 




Fig. 6.— MGRO J2031+41: 1 and 2-a C.L. con- 
tours for the power law fit. The point indicates 
the best fit result (minimum x^)- 



MGRO J2031-I-41, observed with a significance 
of 7.3ct, is an extended source, with a = 1.1°. 
It shows a central core with a higher significance 
(> 6cr) surrounded by a broader region (approxi- 
mately 3° X 4°) with a lower significance (> 4(t). 
TeV J2032-f 4130 lies within the central core, and 
the Fermi Cocoon 1-a radius {cr — 2°) completely 
includes MGRO J2031-F41 However, the Fermi 
analysis ( Abdo et al. 2011b() shows that the Co- 
coon has an elongated shape, poorly correlating 
with the Milagro contours, and the possible as- 
sociation between the Cocoon and the extended 
Milagro excess remains unclear. 

The spectral fit results do not change if the 
maximum significance position is chosen instead 
of the 2-D Gaussian fit position. 

4. Conclusions 

We present the spectra of the two brightest Mi- 
lagro sources in the Cygnus region using a new 
analysis technique applied to the last 3 years of 
data collected by the Milagro experiment. 

MGRO J2019+37 is observed with a signifi- 
cance over 12(7. Its emission is well fitted by 
a power law with an exponential cutoff {Ec = 



29 
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TeV) and a hard asymptotic spectral index 
{a = 2.0tii). The simple power law hypothe- 
sis is disfavoured at a C.L. > 98%. The TeV ex- 
cess measured by Milagro from MGRO J2019-I-37, 
spatially associated with the Fermi LAT pulsars 
J2018.0-H3626 and J2021.0-h3651, has been con- 
firmed by VERITAS, and it is likely produced 
by nearby unresolved sources. ARGO-YBJ does 
not detect a significant emission from MGRO 
J2019-f-37, but the 90% C.L. upper limits do not 
conflict with the Milagro best fitting model. 

The emission from MGRO J20314-41 (7.3cr sig- 
nificance), is well modeled by a power law with 
a — 3.08^q'\^. Our result, in particular the inte- 
gral fiux in the overlapping energy range, is con- 
sistent with previous measurements by ARGO- 
YBJ and Whipple, but it disagrees with HEGRA 
and MAGIC resuhs for TeV J2032-^4130, most 
likely because of the different PSF of the instru- 
ments and different background subtraction meth- 
ods. MGRO J2031-I-41 appears to show a struc- 
tured morphology, produced by the superimposi- 
tion of a central point-like source, coincident with 
TeV J2032+4130, and an extended emission, pos- 
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Fig. 9. — PSF-smoothed significance maps of the Cygnus region, top panel - Black triangles mark the 
positions of the two Milagro sources with the l-cr radius (solid circles), white-filled circles mark Fermi BSL 
sources (the 1-a radius is consistent with the marker size). The Fermi Cocoon position is marked by the 
cross and the dashed ellipse is its l-tr radius, bottom left panel - zoomed MGRO J2019+37 significance map. 
The circle represents the 1-a radius, bottom right panel - zoomed MGRO J2031+41 significance map. The 
Milagro position with the 1-a radius is shown in black. The crosses mark the Whipple (green), HEGRA 
(magenta) and MAGIC (blue) positions with their uncertainty. 




Fig. 10. — Fermi LAT counts maps of the Cygnus region above 10 G eV. We used Pass 7 photon data from 
2008-08-04 15:43:37 to 2011-12-06 10:10:13 GMT (|LAT Data Serveit) . zenith angle < 107°, data_QUAL=1, 
LAT_CONFiG=l. A Gaussian smoothing is applied for display purposes. Circles mark the position of the 
Fermi BSL sources. White lines are the Milagro significance contours (3, 5 and llcr). 
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sibly produced either by unresolved sources or in- 
teractions of cosmic rays with the local interstel- 
lar medium. The correlation between the TeV ex- 
tended emission and the overlapping Fermi Co- 
coon is unclear and needs further studies. 

HAWC, the next generation water Cherenkov 
observatory, will be able to produce a more accu- 
rate analysis of the TeV emission from the Cygnus 
region, with its improved sensitivity (10-15 times 
better than Milagro) and angular resolution. 
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periment. This work has been supported by the 
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Fig. 7.— MGRO J2031+41: 1 and 2-a C.L. con- 
tours for the power law with cutoff fit. left panel 
- NioTeV vs a projection, right panel - Ec vs a 
projection. The points indicate the best fit result 
(minimum x^). The right panel shows that the 
present analysis is unable to constrain the upper 
limit of the energy cutoff within la. 
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Fig. 8.— MGRO J2031+41: energy spectra. The 
power law model is shown in blue and the power 
law with cutoff model is shown in red. These 
two hypotheses give the same fit probabil- 
ity. The pink cross is the Whipple flux at 0.6 



TeV (|Lang et al.l |2004| ). The fine-dashed cyan 



line and the dashed yello w line are the HEGRA 
and MAGIC best fits (lAharonian et al.l l2005t 
2008t ). respectively. The dot-dash 
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black line is the ARGO-YBJ best fit (Bartoh et al 



20121 ). 
band. 



The shadowed area represents the 1-a 
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Table 1 
Source Positions 



survey 


RA 


Dec 


width (ct) 


MGRO J2019+37 








2-D Gaussian with base 


304.646+0.008 


36.847+0.006 


0.75 


Fermi position (OFGL J2020. 8+3649)" 


305.22 


36.83 




Milagro 2007'= 


304.98 


36.66 


1.1 


MGRO J2031+41 








2-D Gaussian witli base 


307.41+0.02 


41.190+0.015 


1.1 


Fermi position (OFGL 12032.2+4122)" 


308.06 


41.38 




Milagro 2007*= 


307.98 


41.51 


3.0 



" lAbdo et al.|[2009h 
'' llAbdo et al.ll2007ah 

Note.— Positions of MGRO J2019+37 (top) and MGRO J2031+41 (bottom) for the 
2-D Gaussian with base fit, Fermi position of associated sources and previous Milagro 
survey. 



Table 2 
Fit Results 



Source 


A^ioToV (10"^" s"^ m"^ ToV"^) 


a 


Be (TeV) 


xVd-o.f. 


prob. 


MGRO J2019+37 
MGRO ,12019+37 


4.1 + 0.6 

' —2 


2.78 + 0.10 


oo 

,0 + 50 


16.12/7 
1.924/6 


0.024 
0.93 


MGRO J2031+41 
MGRO J2031+41 


q 7+90.0 
■'^ ' -1.9 


3.08^S:1? 

9 7+0.5 
^- ' -2.6 


oo 
44+°° 


6.980/7 
5.796/6 


0.43 
0.45 



Note.— Best fit results for MGROJ 2019+37 and MGRO J2031+41 using both the power law and the 
power law with cutoff hypotheses. The table lists statistical errors only, the systematic flux error is ~30% 
and the systematic error on the spectral index is ~0.1. 
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